ABSTRACT OF THE DISCLOSURE 

Domain characterization generated by Voronoi tessellation, which is very close to 
realistic geology and computation of gravity response of such domain, has a three dimensional 
fractal basin structure, and is favorable for oil exploration. Interfaces or tessellating domains are 
represented by a set of parameters, which are referred as Voronoi centers. These parameters can 
be perturbed by any amount without getting into representational problems. To accomplish such 
representation Voronoi tessellation is used, which in two dimensional space involves enclosing 
every Voronoi center by a Voronoi polygon such that the common edge of adjacent polygons is a 
perpendicular bisector to the line joining the Voronoi centers on both the sides of that edge. 
Instead of using conventional Euclidian distances, the notion of Voronoi tessellation is 
generalized by using L p distances, where p can hold any real value so that Voronoi domains are 
not necessarily polygonal. A desired fractal subsurface is generated using this approach that is 
quite close to the natural settings. Next, the gravity response due to this fractal subsurface 
structure is computed. A significant advantage is provided especially in geophysical inversion 
where initial model parameters are updated in each iteration, which can be done more easily and 
efficiently by Voronoi tessellation merely by changing Voronoi centers. 


